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We investigate theoretically the quantum phase transition (QPT) between the one-channel Kondo 
(1CK) and two-channel Kondo (2CK) fixed points in a quantum dot coupled to helical edge states of 
interacting 2D topological insulators (2DTI) with Luttinger parameter < K < 1. The model has 
been studied in Ref.— , and was mapped onto an anisotropic two-channel Kondo model via bosoniza- 
tion. For K < 1, the strong coupling 2CK fixed point was argued to be stable for infinitesimally 
weak tunnelings between dot and the 2DTI based on a simple scaling dimensional analysis^. We re- 
examine this model beyond the bare scaling dimension analysis via a 1-loop renormalization group 
(RG) approach combined with bosonization and re-fermionization techniques near weak-coupling 
and strong-coupling (2CK) fixed points. We find for K — > 1~ that the 2CK fixed point can be 
unstable towards the 1CK fixed point and the system may undergo a quantum phase transition be- 
tween 1CK and 2CK fixed points. The QPT in our model comes as a result of the combined Kondo 
and the helical Luttinger physics in 2DTI, and it serves as the first example of the 1CK-2CK QPT 
that is accessible by the controlled RG approach. We extract quantum critical and crossover behav- 
iors from various thermodynamical quantities near the transition. Our results are robust against 
particle-hole asymmetry for 1 < K < 1. 

PACS numbers: 72.15.Qm, 7.23.-b, 03.65.Yz 



I. INTRODUCTION. 

Quantum phase transitions (QPTs) 1 , the continuous 
phase transitions at zero temperature due to competing 
quantum ground states or quantum fluctuations, in cor- 
related electron systems are of great fundamental im- 
portance and have been intensively studied over the 
past decades. Very recently, nano-systems (in particu- 
lar quantum dots 2 -) offer an excellent playground to study 
QPTs due to high tunability 3 . - — . The well-known Kondo 
effec t 11 ' 12 plays a crucial role in understanding low en- 
ergy properties in quantum dot devices. Potential new 
QPTs in these systems may be realized in connection to 
exotic Kondo ground states. An outstanding example of 
an exotic Kondo state is the two-channel Kondo (2CK) 
syste m 13 ' 14 , which has attracted much attention both 
theoretically and experimentally as it shows non-Fermi 
liquid behaviors at low temperatures. Experimentally, 
the 2CK behaviors have been realized in RefJ^ where a 
quantum dot independently couples to an infinite and a 
finite reservoirs of non-interacting conduction electrons. 

More interestingly, the 2CK physics has also been 
found theoretically in Kondo quantum dot coupled to two 
strongly interacting Luttinger liquid leads with Luttinger 
parameter K < 1 16 i 17 . In this case, electron-electron 
interactions in the leads strongly suppress the cross- 
channel Kondo correlations responsible for charge trans- 
port through the quantum dot while the Kondo correla- 
tions involving electrons on the same lead are unaffected, 
leading to an insulating two-channel Kondo ground state 
where two independent Kondo screenings occur between 
the spins on the dot and in each lead separately. On the 
other hand, for weaker electron interactions, K > i, both 
kinds of the Kondo scattering involving conduction elec- 



trons on the same and different leads become relevant at 
low energies, giving rise to conducting 1CK ground state 
where only a single-channel electrons (even combination 
of the two leads) effectively couple to the Kondo dot. An 
exotic quantum phase transition between the conducting 
1CK phase for K > ^ and the insulating 2CK behaviors 
for K < | is therefore expected at K = However, 
the critical properties of this 1CK-2CK QPT have not 
yet been addressed yet since it is not accessible to any 
controlled theoretical approaches. 

On the other hand, recently a new type of materials- 
topological insulators (Tls)-with a gaped bulk and gap- 
less edge states has been proposed theoretically^ and 
realized experimentally^. In 2D TIs, the gapless edge 
states have "helical" nature, i.e. the directions of 
spin and momentum are locked together—. Based on 
bosonization and a simple scaling dimension analysis, the 
2CK behaviors were argued to be stabilized in Kondo 
quantum dot coupled to two interacting helical edge 
states of 2D TIs as long as a weak electron-electron in- 
teraction exists in the helical electrons (K < l) 21 . How- 
ever, on a general ground, similar competition between 
the cross-channel Kondo correlation and suppression of 
tunneling due to electron-electron interactions mentioned 
above is also expected here for the helical Luttinger liq- 
uid, a special type of Luttinger liquid with broken SU(2) 
symmetry. The exotic 1CK-2CK QPT may therefore oc- 
cur in this new setup. 

In this paper, we re-examine the system in Refill near 
2CK fixed point to explore the possibility of the exotic 
1CK-2CK QPT via the controlled 1-loop renormaliza- 
tion group (RG) approach combined with bosonization, 
which goes beyond the bare scaling dimension analysis 
in Refill. For a weak but finite lead-dot tunneling, we 
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find that for X 4 1~ the 2CK fixed point can be un- 
stable towards the (anisotropic) 1CK fixed point and the 
system may undergo a quantum phase transition (QPT) 
between 1CK and 2CK fixed points. The QPT in our 
model comes as a result of the combined Kondo and the 
helical Luttinger physics in 2DTIs. It serves as the first 
example of the 1CK-2CK QPT that is accessible by the 
controlled perturbative RG approach. The stability anal- 
ysis on these two fixed points shows that they are stable 
against small particle- hole asymmetry for | < K < 1. 
We extract the non-Fermi liquid behaviors from various 
thermodynamical quantities at the 1CK-2CK quantum 
critical point. 

This paper is organized as follows. In Sec. II, we in- 
troduce the model Hamiltonian and its bosonized form 
as shown in Rcf. 21 . In Sec. III. A., we present the 
RG analysis of the model both in the weak coupling 
limit via bosonization approach (see Appendix A.). In 
Sec. III. B., we further map our bosonized model onto 
an effective Kondo model via re-fermionization near the 
strong-coupling 2CK fixed point. We then perform the 
RG analysis via both poor-man's scaling (see Appendix 
B.) and field-theoretical e-expansion technique (see Ap- 
pendix C). Our RG analysis in both limits suggests a 
quantum phase transition between 1CK and 2CK fixed 
points. In Sec. IV. we perform stability analysis on the 
1CK and 2CK fixed points respectively. We find that 
both fixed points are stable for \ < K < 1, which sub- 
stantiates our main finding that there exists an unstable 
quantum critical points separating two stable 1CK and 
2CK fixed points near K = 1. In Sec. V., we calculate via 
field-theoretical e-expansion approach the critical proper- 
ties and crossover functions of various thermodynamical 
observables. In Sec. VI., we emphasize the clear physical 
picture of our main findings and draw conclusions. 



and the electron-electron interactions with forward scat- 
tering g2A > terms are given by Hint with i = 1, 2 the 
lead index, and a — R,L being the label of the right 
(R) and left (L) moving electrons in the helical edge 
state. The conduction electron spin operator is given by: 
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The local impurity spin operator on the quantum dot can 
be expressed in terms of pseudo-fermion operator frr^: 
S = s f\ ' ~2~ ' /*• I n the Kondo limit of our in- 
terest, the impurity (quantum dot) is singly-occupied: 
Eo-=t 1 /a/"- = 1- Here, the couping J\ and J2 in Hk 
stand for the strength of the Kondo correlations between 
the dot and electrons on the same and different leads, 
respectively. Note that in the presence of spin-orbit cou- 
pling, the spins f/4 of the helical edge state electrons are 
locked with their right-moving (i?)/lcft-moving (L) mo- 
mentum. Note also that the small spin-orbit coupling will 
break the SU(2) spin-rotational symmetry in the above 
isotropic Kondo model, leading to the anisotropic Kondo 
model with Jf v ^ J?2L. 

The Hamiltonian Eq. [T] can be bosonized through 
the standard Abelian bosonization 2 ^ for the electron 
operator. 2 ^: c lM/L = -±=F^ R/L e ±l{ -^'>'^'^ xS > +kF ^\ 
the bosonic fields 4>i{x) — 4>iL(x) + <fti R {x), Oi{x) — 
4>iL{x) — 4>in(x). The dual fields (j)i{x) and 9i(x) obey the 
commutation relations: [4>i(x), Oj(x')] = : ^SijSgn(x — x') 
with sgn(x = 0) = 0. The symmetric and antisymmetric 
combinations of <fii, Oi are defined as: 4> s / a — -^("fii i'fe) 
and 9 s / a = ^75(^1 ± 6*2). Here, -F^/l are the Klein fac- 
tors to preserve the anti-commutation relations between 
fcrmions (electrons) in the bosonsized form, and a is the 
lattice constant (lower bound in length scale); we have 
also dropped the spin indices of the edge state electrons 
due to their helical nature. The bosonized Hamiltonian 
after rescaling the boson fields is given by2i: 



II. MODEL HAMILTONIAN. 

In our set-up, the Kondo Hamiltonian has the same 
form as in Ref/ 2 ^., given by: 

H = Ho + Hk + Hi n t , 

Hq = —wf 
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+ .92 J cf R {x)c{ R (x)c%{x)c{ L (x)dx. (1) 

Here, Hq describes the two conduction electron baths (la- 
beled as lead 1 and lead 2) made of helical edge states in 
2D topological insulators, Hk is the Kondo interaction, 
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with K p = l/K a = K = y i+fi« F +i| OT y 

being the Luttinger parameter, v F = 
Vf^J(1 + %£-irvp) 2 — (22-irvp) 2 . Here, we consider 
repulsive electron-electron interactions ($2 , <?4 > 0) , 




FIG. 1: (Color online) RG flows of various Kondo couplings 
in the weak-coupling regime with fixed bare Kondo couplings: 
jTo' Z = 0.01/io and fj,o = 1. Curves with different colors 
indicate the RG flows with different Luttinger parameters K. 
The solid, dashed, dotted, and dot-dashed lines represent the 
RG flows for , jl, j% v ', and j%, respectively. 
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FIG. 2: (Color online) RG flows of various Kondo couplings 
in the weak-coupling regime with fixed K = 0.8 for various 
bare Kondo couplings (in units of /Uo = 1). The solid, dashed, 
dotted, and dot-dashed lines represent the RG flows for j* y , 
jl, 32 V , and j$, respectively. 




giving < K < 1; and K p r a ) refers to the interac- 
tion strength in the charge (spin) sector. Note that 
K a = 1 , K p < 1 corresponds to a spinful Luttinger 
liquid with SU(2) spin symmetry; while as K a ^ 1 when 
this symmetry is broken. The helical Luttinger leads 
in Hq + Hi n t we consider here corresponds to a spinful 
Luttinger liquid lead with broken SU(2) symmetry and 
K a > 1— . Note that we have dropped the Klein factors 
in Eq. [5] as they can be included straightforwardly in the 
same manner as shown in Refsi 22 ' 29 . 



III. 



RG ANALYSIS OF THE MODEL. 



A. RG analysis in weak coupling fixed point: 



0. 



1. RG scaling equations. 

In the vicinity of the fixed point Jf y ' z = 0, it has 
been shown in Refill that the scaling dimensions of these 
Kondo couplings based on the bosonized Hamiltonian Eq. 
1 of Ref.21 are: [J^] = K < 1, [Jf ] = 1, [J* v ] = [Jf] = 
h (K + 4f ) > 1 where K < 1 for repulsive electron inter- 
actions. The authors of Ref. 21 argued that for K < 1, 
under renormalization group (RG) transformations, the 
2CK fixed point is reached as the relevant J\ couplings 
flow to large values with decreasing temperatures; while 
the irrelevant J 2 couplings decrease to zero. However, 
for K being slightly less than 1, K 1~, via 1-loop RG, 
beyond the bare scaling dimension analysis, we find it is 
possible that all four Kondo couplings can flow to large 
values, depending on the values of K and the values of 
the bare Kondo couplings. This seems to suggest the pos- 
sible existence of the 1CK fixed point in the parameter 



space of the model. 

Following Refs i 22 i 23 via the renormalization group 
analysis of the bosonized Kondo model Eq. [21 the 1-loop 
RG scaling equations in the limit of K — > l~ read (see 
Appendix A.): 
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where /1 is the running cutoff energy scale, and the 
dimensionless Kondo couplings are defined as: = 



_ ,.K-1 jxy - z 
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P0J1, 3 2 = PoP 2{ k) J2 ) and 
pofJ.?( K+1 <^~ 1 J$ with [i being a cutoff energy scale, 



and po = -^r = being the constant density of states 
for non-interacting leads (K = 1) and po = 1 being the 
band width of the conduction electrons in the leads. Note 
that the linear term in the above RG scaling equations 
comes from the non-trivial scaling dimensions of the cor- 
responding Kondo couplings; while the quadratic terms 
in Kondo couplings are the corrections at 1-loop order. 

For K — > 1~, both J^' 7 ' terms are marginally irrele- 
vant, \J2 V ' Z \ 1 ~- Therefore, within the validity of per- 
turbative RG, both J^' 2 terms can still flow to a large 
value if bare Kondo couplings J^ v,z are large enough (but 
they are still small, J* v ' z — 0(1 — K) <C 1) or the elec- 
tron interactions in the leads are weak enough, leading 
to (possibly) 1CK fixed point (ie., the quadratic terms 
overcome the linear term in RG equations) . However, for 



small enough bare Kondo couplings (or strong enough in- 
teractions in the leads, K -C 1), J 2 v ' z terms are irrelevant 
and hence the system moves towards the 2CK fixed point. 
As shown in Fig. [TJ in the relatively higher temperature 
(energy) regime 10~ 3 < fi/ fio < 1, with decreasing K 
the system tends to flow to 2CK fixed point where J^ v,z 
flow to large values while J 2 V ' Z decreases with decreas- 
ing temperature (energy). On the other hand, for weak 
enough interactions in the leads, — ^ 1 — , all four Kondo 
couplings tend to flow to 1CK fixed point with large val- 
ues (see Fig. [T]). Similar trend is found for a fixed K — > 1 
and different bare Kondo couplings as shown in Fig. [5] 
It is therefore reasonable to expect a 1CK-2CK quan- 
tum phase transition in the parameter space of Ji V 2 ' z , K. 
However, the weak coupling RG analysis is valid only 
at relatively higher energies, and it breaks down as the 
system gets closer to the ground state, which explains 
the rapid increase of J 2 v ' z in Fig. [T] and Fig. [5] at lower 
temperatures where J± v already exceeds the perturbative 
regime, J± v > 1. In fact, the low energy behaviors are 
determined by the physics in the strong coupling regime. 
Therefore, to address the possible quantum phase transi- 
tion between 1CK and 2CK fixed points, it is necessary to 
be able to access the neighborhood of the strong-coupling 
2CK fixed point as we shall discuss below. 




1-K 

FIG. 3: (Color online) The 2CK Kondo temperature T]f K 
(in units of T%, the Kondo temperature of the corresponding 
non-interacting leads, K = 1) as a function of 1 — K for 
various bare Kondo couplings (in units of fio = 1) via the 
weak-coupling RG analysis. 

B. RG analysis near strong coupling (2CK) fixed 
point. 

1. RG scaling equations and the phase (RG flow) 
diagram 



2. 2-channel Kondo temperature T]f K . 



To probe the crossover between 1CK and 2CK fixed 
points, it is instructive to investigate how the Kondo tem- 
perature Tk changes with increasing electron-electron in- 
teraction in the leads (or with decreasing the value of K 
from 1). Since Jf y becomes more relevant with decreas- 
ing K in the weak-coupling regime ([Jf a ] = K < I), it 
is expected that under RG the system first flows very 
quickly to the vicinity of 2CK fixed point. As shown 
in Fig. |31 we find the Kondo temperature T%P K asso- 
ciated with the 2CK fixed point, defined as the energy 
scale n = Tf K under RG where Jf , Jf « 0(1), in- 
creases rapidly with increasing electron interactions in 
the leads, and its value is much larger than the Kondo 
temperature of the same setup in the non-interacting 
limit (K = 1) 1%, Tff K > Tg-. By contrast, in the 
case of a Kondo dot coupled to ordinary spinful Luttinger 
liquid leads in Refsi 16 ' 17 , Jf y is a marginal operator at 
tree-level ([Jf] = 1) in the weak-coupling limit; there- 
fore, the 2CK energy scale T%F K is much smaller than 
the Kondo scale for the corresponding non-interacting 
leads T£, TjP K < Though the system in the weak 
coupling regime quickly approaches the strong-coupling 
2CK fixed point as fj, — > Tjf K , the ultimate fate of the 
ground state depends on the RG flows of various Kondo 
couplings in the strong coupling regime as discussed be- 
low. 



The authors in Refj^i performed scaling dimension 
analysis near a strong coupling regime where Jf = 
0(1),J 2 V ' Z = Jf = 0. They performed the Emery- 
Kivelson unitary transformation^ JJ = e «v27r0<, ^ ^ e 



bosonized Hamiltonian (Eq. 
Eq. 2 of Ref,2i. 

H = H -j?S z d x 9 s (0) 
2 J? 



1 of Refill), and arrived 



^S z sin(W — 9 a (0)) sin(V2^iG/>a(0)) 
ira V K 



(S- + S+)[-±- cos(V2^> a (0)) 
ira V K 



(4) 



with -Jf = V2nKv F -y/^Jf- 

They found the scaling dimensions for the Kondo cou- 
phngs to be [Jf] = f , [Jf ] = ^ [J|] = i {K + i). 
The Jf term is relevant for K < 1, and Jf term is rele- 
vant for h < K < 1 . As K — > 1 ~~ , Jf becomes marginally 
irrelevant, and it can flow to a large value under RG if the 
bare Kondo couplings are large enough once the 1-loop 
RG is performed. This seems to suggest a stable 1CK 
near strong coupling regime as all of the four Kondo cou- 
plings can either flow to or stay at large values (of order 

I)- 

To gain more insight into the stability of the 1CK/2CK 
fixed point, we apply RG approach at 1-loop order 
together with bosonization and re-fermionization near 
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2CK fixed point. First, we shall map the bosonized 
Hamiltonian Eq. 1 of Rcf. 21 onto an effective Kondo 
model via re-fermionization. It has been shown in Ref. 21 
that near the strong coupling 2CK fixed point J xy — > 



oo, Jf -► £>(1),J* 
reads: H 2 ck = H - 



/xy 



► 0, the effective Hamiltonian 
S x cos(V2^K(j) a (0)). Note that 



near 2CK fixed point, the dominating "backscattering" 
J^ y term effectively cuts the Luttinger wire into two sepa- 
rate pieces at x = 21 ' 23 , leading to the well-known open 
boundary condition for an impurity in a Luttinger liq- 
uid at x = 0: c i;R (0) = -c itL (0) (or 4> iR {Q) = -<j) lL (0)). 
The boson field <fi a is approximately pinned to a constant 
value^i. Also, since S x commutes with H2CK, we may 
therefore set S x to its eigenvalue ±5 in Hick- The the 
scaling dimensions of Kondo couplings near 2CK fixed 
point are2i^> [J*v] = i, [Jf ] = 1 + [J|] = i + f . 
Note that all the above three couplings are irrelevant for 
K < 1. This suggests that the system favors the 2CK 
fixed point at ground state K < 1. Meanwhile, by a sta- 
bility analysis in Sec. V., we show that the 2CK fixed 
point is a also a stable fixed point for K > \ once the 
system gets there. 

However, as suggested in our weak-coupling RG anal- 
ysis, the 2CK fixed point may be unstable for K — >• 1~ 
and/or large enough bare Kondo couplings such that 
J^P' 2 may become relevant again, and the system can 
undergo a 1CK-2CK quantum phase transition. To ad- 
dress this possibility, we shall focus below on the 1-loop 
RG flows of the leading two irrelevant operators near the 
2CK fixed point, given by: 



SH 



j*y HTZ 

2CK = -^S x cosU — O (O)) - J?SA9 S (0). (5) 
wa V K 



Via the similar re-fermionization as shown in Appendix 
A., we map Hq + SH2CK onto an effective Kondo model 
subject to a bosonic environment: 



where the boson field Q' a in Hq is decoupled from Hq 
and is added here just for the mapping, the effective 
non-interacting electron operator c% i is defined in Eq. |9] 

Note that since the scaling dimension of cos(^f^-6 a ) at 

2CK fixed point in Eq. |6] is -4 due to open boundary 
condition 2 ^, we have made the following decomposition 

for the boson field \J j^0 a : 



8s = #0, s, 



K 



-w a . 



(7) 



The re-fermionization of Hq is done through the following 
identifications: 

V2# ,a = V2{4 tl + 4 t2 ) = -V2(4 }2 + 4 tl % 

V29 0iS 



4,1 -00,1 + 00,2 ~ 00,2 



(8) 



where we have decomposed the boson field y 8 a into 
two independent sets of boson fields: the "free" (#o,a) 
and "interacting" (9 a ) parts: Here, the "free" part of the 
boson fields $o,a (defined in the same way as in Sec. II.) 
can be re-fermionized into two effective non-interacting 
fermion leads described by Hq with cf a being the electron 
destruction operator of the effective non-interacting leads 



gtU) 
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Hq + SH 2 CK — > Hq + Hq + 8H2CK 
= Hq + H 6 + H 2 CK, 
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dx{d x e' a 



1 \2 



-f J dx[2(d x 6 , a ) 2 + (d x 9 s ) 2 
(d x <p s ) 2 + (d x cb a ) 2 } 

k,cr,i=l(L),2(R) 



H b = ^ j dx2(dj a ) 2 , 
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-Ji(s z LL + s z M )S z 
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with a = l(L),2(i?) being the index for effective non- 

l~fQ^±(z)~5 



interacting leads, s^jf^ 



-k'/3 



being 



the spin- flip (z-component of the spin) operators between 
the effective leads 7 and (3. Note that the effective non- 
interacting leads also exhibit the helical nature; namely, 
the spin up/down (er =f / l) electrons are tied to the 
right (R)/left (L) moving particles, respectively. The 
"free" part of boson field 6>o, Q follow the correlations of 
the free fremions in ID: 



1 

> oc -. 

t 



(10) 



Meanwhile, Hi, represents for the effective dissipative 
ohmic boson environment (baths) made of the "inter- 
acting" part of the effective bosons 6 a . These bosons 
couple to the Kondo dot through the additional exponen- 
tial "phase" factors in the effective Kondo terms H2CK, 
leading to all the combined Kondo-Luttinger physics 29 . 



J 1 

2CK 

FIG. 4: (Color online) Schematic diagram of the RG flow 
near 2CK fixed point. The 1CK-2CK quantum critical point 
(QCP) is represented by the filled black square located at 



In particular, since these dissipative ohmic bosons obey 
the following correlations via Eq. O 



< e 



-iy/ln(±-l)§ a (t) ^4^(^-1)^(0) 



>cx 



1 



(11) 



while the impurity spin operator S z exhibits the following 
correlational: 



<S z (0)S z {t) >oc^. 



(12) 



These correlations lead to the non-trivial bare scaling di- 
mensions of the Kondo couplings and therefore to the 
first term (linear in the Kondo coupling) of the RG scal- 
ing equations. Note that since near 2CK fixed point 
4> St a fields are decoupled from Eq. O we have effectively 
two independent degrees of freedom left among the four: 
{4> s .a,0s,a,)] the open boundary condition for the spin-up 
right-moving (R) and spin-down left-moving electrons: 
cg^ = -cj)^f } is implied in Eq. E With the help of 

Eq. [7] and Eq. [FJ , we finally arrive Hq and H 2 ck in 
Eq.© 

Next, we shall obtain the one-loop RG scaling equa- 
tions for J% v and Jf in Eq|51 To this aim, we define 
the dimensionless couplings j% v = pQC 2 pf J% v (p) and 
jf = poc\pf ' Jf (p) and e = i - 1, e' = 5^ with cf, 
c 2 being defined in Appendix B. and C. 

We derive the 1-loop RG scaling equations via the 
poor-man's scaling approach in Refj^ (see Appendix B.) 
and via field-theoretical e-expansion technique (see Ap- 
pendix C): 
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d In// 



■xy -xy -z 

= m ~ 3i 3i, 



(13) 



For K — > 1 _ we find an intermediate quantum critical 
fixed point (QCP) at j c = (ji c ,j%%) — (e, \ZeP) within 
the validity of the perturbative RG separating the 1CK 
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FIG. 5: (Color online) RG flows of various Kondo couplings 
in the strong-coupling regime near 2CK fixed point with the 
following parameters: (a). j?$*(jx = T]f K ) = 0.4/i , K = 0.8 

K 
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and (b). j«g. 



0* 



0.1/tto, 



0.7. The RG flows 



starting from p = T K to p — > in (a) are towards the 
1CK fixed point; while as the flows in (b) are moving towards 
the 2CK fixed point. The solid, dashed, and dot-dashed lines 
represent the RG flows for jf , j% v , and j£, respectively. Here, 
we set fio — 1. 



fixed point (for j = (ji,o>32fi) > 3c = UtcJZ)) and 
2CK fixed point (for jo < j c ) where j% y and jf flow to- 
wards a large and vanishingly small value, respectively 
(see Fig. S|). The RG flows near the QCP are deter- 
mined by linearizing the RG scaling equations as shown 
in Fig. |4j The typical RG flows corresponding to the 
1CK and 2CK fixed points are shown in Fig. [5] (a) and 
(b), respectively. 

Note that jf — > near 2CK fixed point suggests that 
the original Kondo coupling Jf (see Eq. [I]) is at a large 
value (order ofl): Jf « 0(1), consistent with the famil- 
iar 2CK fixed point with both J* v ' z being large. How- 
ever, we find the "1CK" fixed point here in the strong- 
coupling analysis seems somewhat different from the fa- 
miliar (conventional) 1CK fixed point we obtained in 
the weak coupling regime where all the four Kondo cou- 
plings will flow to (or stay at) large values. Instead, our 
RG analysis based on re-fermionization for the coupling 
jf = Poc 2 p 1 < + 1 <~ 1 J% m Eq.|4]near 2CK fixed point shows 
that it stays irrelevant up to 1-loop order with the RG 
scaling equation: 



djj 
(9 In /1 



1 K 



32 



(14) 



where we find no corrections at 1-loop order. Note that 
unlike in the weak coupling RG where ji V j 2 v will con- 
tribute to the 1-loop renormalization of jf (see Eq[3l the 
ji V j 2 y term is absent here in Eq.[T3]as is already very 
large near the 2CK fixed point, j* v (T w T]f K ) > 1. 

It is clear from Eq. [TJ] that^i j^ip) oc van- 
ishing as p <C Tjf K even for jo > j c where the system 



eventually flows to the ICK fixed point (see Fig. [5]). By 
combining the 1-loop RG analysis in the weak and strong 
coupling limits, we may obtain the full crossover of Jf for 
the system which will eventually flow to the fCK fixed 
point: For T%P K < T < fi , jf first grows to order of 
f (see Fig|TJ); then it vanishes in a power-law fashion at 
lower temperatures T -C T^ CK (see Fig|5]). However, the 
above qualitative feature for jf based on the 1-loop RG 
analysis might get modified at the 2-loop order, which ex- 
ceeds the scope of our current work and will be addressed 
elsewhere. 

Though somewhat unconventional, the fixed point 
with ii2)if — > o° an d jf — > can still be regarded 
as the one-channel Kondo (ICK) fixed point since the 
two leads are connected by the strong transverse Kondo 
couplings j 2 , and only one channel of conduction elec- 
trons (even combination of the two leads) couples to the 
Kondo dot. Therefore, we expect the linear conductance 
G±(T) contributed from J% y at the ICK fixed point here 
to show the same temperature dependence as those in 
the isotropic one-channel Kondo system. 



2. 1-channel Kondo temperature T K 



ICK 



As mentioned above, for jo > j c with decreasing tem- 
perature the system crosses over from 2CK to ICK fixed 
point at a much lower energy scale fi ~ T]f K <C Tff K 
where T]^ K refers to the Kondo temperature associated 
with the ICK fixed point. As shown in Fig. [6l the ICK 
fixed point persists to be the ground state at a finite but 
weak electron-electron interaction strength, K c < K < 1 
with K c being the critical interaction below which the 
ground state switches from ICK to 2CK fixed point. 
Meanwhile, the crossover scale to ICK fixed point T]f K 
(with respect to T]f K ) for jo > j c gets reduced sig- 
nificantly as interaction gets stronger. Also, for a fixed 
value of K, the ratio T^P K /Tjf is larger for larger bare 
Kondo couplings jo, as expected. 



IV. STABILITY ANALYSIS OF ICK AND 2CK 
FIXED POINT FOR K < 1. 

Having found the possible QPT between ICK and 2CK 
fixed points, it is important to perform a stability anal- 
ysis and study how robust the quantum critical point of 
our system is against small perturbations. Equivalently, 
we need to know how stable the ICK and 2CK fixed 
points are for K < 1. 

We first examine the stability of the helical Luttinger 
liquid lead itself. In general there exists the single parti- 
cle backscattering term due to the interaction of c iR /^(0) 

and the quantum dot 21 : t'c\ R Ci t L + h.c. However, this 
term is forbidden here as it breaks time-reversal sym- 
metry. Meanwhile, for 1-D Hubbard model in general 
there exists the "spin-flip" backscattering term in H u of 
the form: H s f oc c^ L c\ L c\^ R c\ R + h.c. However, due to 



r r2CK 
K 0.4 




FIG. 6: (Color online) The ratio of ICK to 2CK Kondo 
temperature Tjf K /T]f K as a function of 1 — K for various 
bare Kondo couplings j^o' z (/i = fJ-o) (in units of fio = 1) via 
the weak-coupling RG analysis. 



the helical nature of our leads (or the Right /Left moving 
electrons are tied to their spins, i.e., only c]yn electrons 
exist), this H s f term is therefore absent. Nevertheless, 
the the Umklapp term that exists on a single bond is 
allowed by the time-reversal symmetry 2 ^: 



= ^^(0)^(0)4(0) 

x 4(0) + h.c. 



(15) 



The scaling dimension of this term has been shown to 
be [H um ] = 4K, suggesting that the helical edge state is 
unstable towards an insulating phase for K < j. 

We now focus on the effects of the particle-hole (p- 
h) asymmetry on the stability of these two fixed points 
as indicated in Refs i 16 i 17 ' 29 that it is the most relevant 
perturbation for a Kondo dot coupled to Luttinger liquid 
leads. Let us first address this issue at 2CK fixed point 
where the two leads are effectively disconnected. The 
particle-hole asymmetry in our Kondo model generates 
potential scattering terms of the following formi^: 



H 



IpS 



Ht + H te , 
t 



t, 



E J<?a act 

k,i,a,a=L,R 

EJaa aa 
C k,i c k,j 



h.c. 



(16) 



with i,j = \{L),2(R) being the lead index, a =f (R),i 
(L) being spin index, and R(L) being the right (left) 
moving particles. Here, t and i e terms represent a chem- 
ical potential of each lead and a weak tunneling between 
the disconnected Luttinger leads 2 ^. Meanwhile, two ad- 
ditional two-particle scattering terms H2 P involving tun- 
neling of spin (to-) and of charge (t p ) can be generated 
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by the weak tunneling t e via 2nd-order perturbation (see 
Fig.2 (d) (e) (f) of Ref,2S), given by: 



H- 



2ps 



H, 



H t +H t , 

L ° 2-^i M L k,2 L k,2 L k,l ' ll - c - 



C k,2 C k,l C k,2- 



The bosonized form of Eq. [19] reads 28 : 
K 



Hip S + H; 



2ps 



(17) 

(18) 
(19) 



t- — d x d> s 
2na Y 

cos(V2^ o (0)) cos(W ^0„(O)) 
Z7ra V K 



and i> p terms are irrelevant for ^ < K < 2 and relevant 
otherwise. 

Based on the above analysis, we find that both 1CK 
and 2CK fixed point are stable for h < K < 1, and un- 
stable for if < |. We have checked that our analysis 
reproduces the well-known results for a Kondo dot cou- 
pled to conventional Luttinger liquid leads in Refs i 16 ' 17 i 29 
where [t e ] = at 2CK fixed point and [v e ] = i(l + K) 
at 1CK fixed point. 

As a final remark, we consider here the parity (left- 
right) symmetric model where Ji = Jll = Jrr with 
Jll(rr) being referred to the Kondo couplings involving 
only the left (right) lead. Nevertheless, parity asymme- 
try is a relevant perturbation near 2CK fixed point. In 
the presence of parity asymmetry (Jll ^ Jrr), the sys- 
tem will flow to the 1CK fixed point with the large bare 
Kondo couplings 2 ^. 



+ cos(2V2^</) a ) 
2na 



(20) 



Near 2CK, (f> a (0) is a constant, therefore the scaling di- 
mensions of these term gives: [t] = l,[t e ] = t^, [ta] = 
(t p w const.)2&. It is clear that all operators are irrele- 
vant for | < K < 2; t e term becomes relevant for K < i, 
and t a is relevant for K > 2. 

Next, we consider the stability of the 1CK fixed point. 
Since the cross-channel Kondo coupling J% v term flows 
under RG along with J* v to large values while as Jf 
stay at order of 1, the two semi- infinite Luttinger wires 
are joined into one single infinite Luttinger wire 7 . In 
contrast to the "weak tunneling" processes mentioned 
above at the 2CK fixed point, the potential scattering 
term generates the "weak backscattering" processes be- 
tween the electrons in the upper and lower edges, includ- 
ing the single-particle backscattering term v e , and the 
two-particle backscattering terms v p , and v„ (see Fig. 2 
(a) (b) (c) in Re&2£): 



H v 



v V c n V L + c UL c iR + h c 



= v P y:^^^+h.c, (2D 



k 

V v C k,l C k,2 C k,2 C k,l ' n - C - 



(22) 
(23) 



In fact, there exists a duality mapping between the "weak 
tunneling" and "weak backscattering" limits 



JR 

"k,2 



L fc,2' L fc,2 



-> C 



V P .t, 



Note that at 1CK fixed point, <\> a is not pinned to a con- 
stant as opposed to that in the 2CK case. The scaling 
dimensions of these terms can be read off straightfor- 



wardly: [v e ] = \[K + i), [v p ] = 2K, and [v a ] 
The v e term is always irrelevant for K < 1, while the v a 



_2_ 

K ■ 



V. CRITICAL PROPERTIES NEAR 1CK-2CK 
QUANTUM PHASE TRANSITION. 

The critical properties and crossovers of various ther- 
modynamical quantities near this newly found 1CK- 
2CK QCP can be obtained via the above RG ap- 
proach combined with the field-theoretical e— expansion 
technique 3 -^—. We employ here a double- e— expansion 
with two small expansion parameters e and e'. Our ap- 
proach is is valid for the Luttinger parameter K — > 1~ as 
both parameters e and e' are within perturbative regime: 
e -> 0, e' \ < 1, and e < e'. Following Refs.— ~— , 

we define the renormalized pseudo-fermion fields f a and 
the renormalized dimensionless Kondo couplings j as: 



/a 



and J 2 y 



■XV 

-32 . 



with and / 2 being the renormalization factors 
for the impurity field and Kondo couplings, respectively 
and n is a renormalization energy scale. The renormal- 
ization factors are obtained via minimal subtractions of 
pole a 32 i 33 , given by (see Appendix C): 



Zf = 



e'' 

2e ' 

(Jl y ) 2 , (jT) 2 

8e 16e' ' 



(24) 



Within the field-theoretical RG approach, we have 
checked that the RG scaling equations in Eq. Q2] 
can be reproduced via calculating the /3— functions: 
f3(ji) = Ma^rliio w ith fi being an energy scale, 
ji = j^jji being the renormalized Kondo couplings 
and j 1>0 = jf >0 = Jf,j2,o = j'2,0 = J 2 V bein S the bare 
Kondo couplings (see Appendix C). Below we discuss 
various critical properties and crossover functions based 
on field-theoretical e-expansion approach. 
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A. Observables at criticality. 

We first calculate various observables at criticality, 
including correlation length exponent, impurity entropy, 
dynamical properties of the T-matrix and local spin 
susceptibility. 



1. Correlation length exponent v . 

The correlation length exponent v describes how the 
correlation length £ diverges when the system is tuned to 
the transition: £ oc \t\~ v with t = ''"J 3 " being the dimen- 
sionless distance to the QCP. It also gives the power- 
law vanish of the characteristic crossover energy scale 
T* close to the transition: T* oc \t\ u . To calculate v, 
we first linearize the RG scaling equations EqQ~3] near 
QCP. The correlation length exponent v is determined 
by the largest eigenvalue of the coupled linearized equa- 
tions, found to be: 



VI + 16ife-l 2e 



(25) 



where the leading order behavior v sa ^ is obtained by 
expanding the square- root in Eq.[5I)]in the limit of e -C e'. 



2. Impurity entropy. 

The impurity contribution to the low-temperature en- 
tropy near QCP is obtained by a perturbative calculation 
of the impurity thermodynamic potential Q imp ££ with re- 
spect to the 2CK fixed point and taking the temperature 



derivative: SV, 
written as: 



dT 



At QCP and T — it can be 



S 



QCP 
imp 



s 



2CK 



AS, 



imp ■ 



(26) 



In \/2K = \ In 2K is the zero-temperature 



where S?£? 

residual impurity entropy at 2CK fixed point which 
shows the existence of fractionally degenerate ground 
stateM£L22, and AS lmp is the correction to Sf^ at 
QCP. Following similar renormalized perturbative calcu- 
lations in Reb 3 ^, we find 



ASi m p = it 2 ln2[ 



3e 2 7T 2 In 2 
32K 



(27) 



scattering of the conduction electrons from lead a to 
lead a' via the impurity. In particular, T aa i (lo) with 
a 7^ a' describes the transport across the dot, de- 
tectable in transport measurements. The T— matrix is 
determined from the conduction electron Green func- 
tions G aa '(t) =< c a (0)c' a ,(t) > through G aa i = 
G° aa ,S aa , + Gi a {Lo)T aa ,(u)Gl, a ,{u;^.. Near 2CK fixed 
point, T aa '(u>) with a ^= a' is defined through the 
propagator, Gt, of the composite operator T aa = 

j-y e WMi-i)Sa(o) f t u ,^ (see Eq . m . Taa , [u) = 
G T (w)aa. 

Following the similar calculations in Ref<22 and Ap- 
pendix C, we analyze the propagator Gt(oj) near 2CK 
fixed point and find at zero temperature Im(T 2 ^, K (u))) oc 

_ \ck with the anomalous exponent at the tree level 

with respect to the 2CK fixed point given by j] T CK = 
2e (ie., Im(T^, K (w)) oc w 2e ). Near 1CK-2CK QCP, 
however, Im(T aa >(uj)) acquires an additional anomalous 
power-law behavior: 



Im(T aa i(u))) oc 



UJ 



(29) 



where the additional anomalous exponent r/T is obtained 
via the renormalization factor Zt for the T— matrix 
propagator T aa ,{oj^M*: rfr = /3(jf j tc + 

P(Ji ) 9 l (>f* T I oV >il ' Here, the renormalization factor Zt is 

obtained by minimal subtraction of pole s 32 ' 33 : Zt = -^f- 

nr-n ^ 

with Zf, Zj± given by Eq. [Ml We find therefore 



{Jit? 



2& 



(Jfc) 2 



e 

Ik 



2 4 
and Im(T aa i(uj)) at QCP behaves as: 

Im(T aa r(uj)) oc ujTk+— 



2e 



(30) 



(31) 



4- Local spin susceptibility \ zz (uj). 

The local dynamical spin susceptibility lTn(Xzz(&)) a t 
the impurity (quantum dot) is defined as the time Fourier 
transform of the spin-spin correlator: < S z (0)S z (t) >. 
At zero temperature, the imaginary part of the local sus- 
ceptibility, Im(Xzz(w)), shows a power-law behavior at 
QCP: 



Therefore, we have: 



Im(x(uj)QCp) oc — 



1 



(32) 



S 



QCP 



1 nr/ 3e 2 7T 2 ln2 
-\n2K+ . 

2 32K 



3. The T— matrix. 



(28) 



The conduction electron T— matrix, T aa '{uj), in the 
Kondo model carries important information on the 



Here, ifP = e is the anomalous exponent of Im(xzz( u> )) 



at the tree level with respect to the 
point via the correlator < S z (0)S z (t) >oc 



2CK fixed 
-T- evalu- 

tK 

ated at 2CK fixed point^i (ie., Im(x 2 zz K (ui)) oc w £ ), 
and r\ x is the correction to the anomalous exponent 
r] x CK when the system is at QCP. Via e— expansion 
within the field-theoretical RG framewor k 32 ' 33 , r/ x reads: 
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quantum 



critical 



2CK 



T 



1CK 



t = G-i)/j r 



FIG. 7: (Color online) Schematic finite temperature phase 
diagram of our model near 1CK-2CK quantum critical point 
(QCP) located at t = (jo = jc)- As T -> 0, the 1CK 
ground state is reached when t > or jo > j c ; while as the 
2CK ground state is reached for t < or jo < jc- Here, 
jo > (<)jc refers to the upper right (lower left) region in the 
phase diagram shown in Fig.|4]where the RG flows are towards 
the 1CK (2CK) fixed point. The dashed vertical arrow at 
t < refers to the finite temperature crossover between the 
quantum critical region (blue shaded area bounded by the 
crossover temperature T*) and the 2CK ground state, which 
is our interest. Here, jo and j c are defined in the text. 



7 and 7' via the T — > limit of the local susceptibility 
near criticality are defined a o 34 i 35 : 

Xioc(t <0;T = 0) cx (-t)-^, 7 = i/(l-r?x), 
T X ioc(t > 0;T = 0) oc f\ i = vr] x , (36) 

Meanwhile, the critical exponents j3 and 6 associated with 
the local magnetization m; oc can be determined by: 



1 



m loc (t >0;T = 0) cx t p , p = -ur] x , 
m loc (t = 0,T = 0) cx \h\i, S= — -l. (37) 

With the values for critical exponents v (Eq. |2"S"1) and r\ x 
(Eq. [33]) at hand, the other critical exponents are there- 
fore given by: 



7 
/ 

7 
5 



1 l + Ke 



2e 

1 + ife 



+ 0(e 2 ,e' 2 ), 



8K 
l + Ke 
16K 
8K 



Ke 2 



+ 0(e 2 ,e' 2 ), 
+ 0(e 2 ,e' 2 ), 
-l + 0(e 2 ,e 12 ) 



(38) 



with Z x = 



Z'j being the renormalization factor for the impurity 
susceptibilit y 32 ' 33 and Zf defined in Eq. [2U Carrying 
out the above calculations, we arrive at 



Vx 



e 



4 ' 



(33) 



and finally Im(x 2( ^ K (uj)) at QCP shows the following 
power-law behaviors: 



B. Hyperscaling. 



(34) 



The impurity correlations at QCP is expected to obey 
certain hyperscaling properties. For example, the local 
dynamic spin susceptibility at criticality obeys y scaling 
in the following form: 



Im(xioc(u,T)) 



A 



_ ri 2CK 



-Vx T' 



(35) 



with being an universal crossover function for the 

QCP here and A being a non-universal pre-factor. Simi- 
lar scaling form can be found in the T-matrix. Hyperscal- 
ing can be used to determine relations between various 
critical exponents. It has been known^Sr— that the cor- 
relation length exponent v and the anomalous exponent 
r\ x are sufficient to determine all critical exponents asso- 
ciated with a local field h. In particular, the exponents 



C. Crossover near critical point. 

Next, we focus on calculating the crossover functions 
close to the 1CK-2CK quantum critical point. In general, 
the crossover functions of observables near criticality de- 
pend on the RG flows of both jf and j% v (see Fig. [7]); 
therefore they may not be expressed analytically in terms 
of universal crossover functions of a single variable. Nev- 
ertheless, great progress can be made when one makes a 
special choice of bare (initial) values of Kondo couplings 
such that e'ji o = (j% o) 2 - Note that this set of bare cou- 
plings can in general be tuned through adjusting various 
microscopic parameters, such as: spin-orbit coupling in 
2DTIs, the lead-dot hoping. For this particular choice 
of bare couplings, we found and checked that the RG 
flows of and j^Qji) follow the well-approximated 

trajectory: e'jf ~ (j^) 2 > {i-e., fl(jf) ~ 0). Under this 
constraint, only one RG j3— function {^(j^)) effectively 
remains: 



/3(j2 y ) 



■ xy 



2K(j 



xy\3 



(39) 



One can therefore easily solve Eq. [39] analytically, and 
its solution for the range between QCP at = \fS = 



ee' and the 2CK fixed point (j 2 < j 2 c where our RG 
and e-expansion approach is controlled) is found to be: 



■ xy 

32c 



(40) 



where T* = ( 



EST" 



) 2t is the crossover energy 



scale. It is clear that the power-law vanish of T* follows: 



I 




FIG. 8: (Color online) X™ P /x?mf versus T/T* (see Eq|2If 
at fixed J£" = 0.8 for various bare Kondo couplings (in units of 
Ho = 1) where jf = 2K{j2 V ) 2 is satisfied and = ~ 
0.39. 
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FIG. 9: (Color online) Crossover of the linear conductance 
G(T) versus T/T* (see Eqi3} at fixed K = 0.8 for various 
bare Kondo couplings (same as in Fig(8}. 



T* cx |t|A = \t\ v with the correlation length exponent 
v being v — which agrees with our earlier result in 
Eq. [25l The crossover function in Eq. [40] can be used to 
compute various crossovers in thermodynamic functions 
near 1CK-2CK QCP as discussed below. 



The impurity susceptibility Txim P (T) . 



susceptibility is 



,33.34. 



defined as 

^ where 



The impurity 

Ximp (T) = Ximp,imp + ^Xu.imp + (Xu.u ~ X^u) 

Xu,u is the bulk response to the local field applied to 
the bulk only, Ximp,imp is the impurity response to the 
local filed applied to the impurity only, Xu,imp is the 
crossed response of the bulk to an impurity field, x b u d u 
is the susceptibility of the bulk in the absence of the 
impurity. We can calculate Ximp(T) via perturbative 
approaches in Refs. 33,34 . We find (up to the first order 
in j^) Ximp(T) has the following crossover form (see 
Eq. inland Fig.©: 



Ximp 

(T) 

Ximp 



1 - 



■ xy 
32c 



i + (£) 



-2c 



(41) 



where xfmp 1S the impurity susceptibility at the 2CK 



act 



dT 



fixed point, given by* 3 - Ximp {T) oc 

with impurity specific heat at 2CK fixed point given 
by: Cf£« cx Tk- 2 (for § < K < 1) and C%£ cx T 
(for K < |)2L22. We have therefore rf^* = -| - 2 (for 



1<K<1) 



and i] 



2CK 



1 (for § < if < 1) 



2. Impurity entropy Si mp {T). 



At 2CK fixed point, the impurity residual entropy has 
been calculated in Refi2I: Sf£p = \n\/2K. Following 
Ref.— , the correction to Sf^f near QCP is obtained 
within perturbative RG approach the by calculating the 
thermodynamic potential and taking the temperature 
derivative. The crossover function for the impurity en- 
tropy near QCP is found to be^ 3 -: 

Simp (--0 



C2CK 
imp 




3. Equilibrium conductance G(T). 

The equilibrium conductance G(T) has the following 
crossover form between 2CK fixed point and the QCP 
(see Fig. HI): 



G(T)oc [£»(/, ->T) 



T 2e + (T*) 2 



(43) 



Note that in equilibrium the linear conductance at 2CK 
fixed point G2Ck(T) is determined by the bare scaling 
dimension of the leading irrelevant operator j 2 v , \j% ] = 
This gives G 2C k{T) oc T 2 ^- 1 ) = T 2e . For T < 
T* where the system reaches the 2CK fixed point, the 
temperature dependence of G{T) in Eq. 03] reduces to 
that at 2CK, G(T < T*) cx G 2C k(T), as expected. 



VI. DISCUSSIONS AND CONCLUSIONS. 

Before we conclude, we would like to emphasize again 
the clear physical picture we provided in the Introduction 
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to make our main results more transparent. First, it is 
well-known that the stable one-channel and two-channel 
Kondo fixed points are expected in the case of Kondo 
quantum dot coupled to two conventional spinful Lut- 
tinger liquid lead s 16 ' 17 . The ground state of this system 
changes from 1CK to 2CK when Luttinger parameter 
K reduces from the non-interacting limit (K — 1) to the 
strongly interacting limit (K < i). The electron-electron 
interactions in the Luttinger liquids act equivalently as if 
an additional dissipative Ohmic boson bath is coupled to 
the quantum dot&22, leading to suppression in electron 
transport from one lead to the other through the dot. A 
quantum phase transition between 1CK and 2CK fixed 
points was argued to exist at K = | as a direct con- 
sequence of the competition between the cross-channel 
Kondo coupling Jlr and the suppression of tunneling 
due to electron-electron interactio n 16 ! 17 . However, up to 
now there is no analytic and controlled approach to ac- 
cess this transition. Note that the 1-loop RG approach 
does not work here to reach to the 1CK-2CK quantum 
critical point since near 2CK fixed point the Kondo cou- 
plings Jll/rri involving in the renormalization of the 
cross-channel Kondo coupling Jlr, both go to infinity 
under RG. 

When a quantum dot couples to helical Luttinger liq- 
uids (a special type of Luttinger liquid), we expect the 
1CK and 2CK ground states are also the two possible sta- 
ble phases for the same reason mentioned above. How- 
ever, due to the helical nature of the Luttinger liquid 
leads, the underling two-channel Kondo model becomes 
anisotropic (Jf y ^ Jf ) as the SU(2) symmetry of the 
model is broken; while as the Kondo model is isotropic 
(Jf y = Jf) for a quantum dot coupled to conventional 
Luttinger liquid leads. This crucial difference enables us 
to access the QPT between 1CK and 2CK fixed point of 
our system via the controlled RG approach. 

In the limit of a weakly interacting helical liquid K — > 
1~ , we find the similar competition between these two 
possible ground states. The 1CK phase is reached when 
J% v is large enough; while as the 2CK phase is reached 
when the electron-electron interaction becomes strong 
enough. Via a controlled perturbative RG approach at 
1-loop order, we find that the 1CK-2CK quantum phase 
transition occurs near K = 1~ . To reach the 1CK-2CK 
phase transition in our setup, we believe it is necessary 
to go beyond the tree-level bare scaling dimension anal- 
ysis, which predicts a stable 2CK phase for as long as 
K < l 2 i. The 1-loop RG is the leading correction to the 
above-mentioned bare scaling dimension analysis. Note 
that the main difference between the case for conven- 
tional Luttinger liquid and that for helical liquid is that 
the resulting two-channel Kondo model is isotropic in the 
former case; while it is anisotropic in the latter case. This 
difference affects details of the critical properties, such as: 
the critical points occurs at K = K c = i for the Kondo 
dot coupled to Luttinger liquid; while as K c = 1 — 0(e) 
for the case of helical Luttinger liquid. At a general level, 
however, we should expect a 1CK-2CK quantum phase 



transition to exist in both cases. 

In summary, we have re-examined Ref^i on the two- 
channel Kondo physics in the Kondo quantum dot cou- 
pled to two helical edge states of 2-dimensional topo- 
logical insulators. Via the 1-loop renormalization group 
approach which goes beyond the scaling dimension anal- 
ysis in Refi^i, we found the quantum phase transi- 
tion between the one-channel (1CK) and two-channel 
(2CK) Kondo ground states for weakly interacting leads 
(K — > 1~). We made definite predictions on the criti- 
cal properties when the system is close to the transition. 
Our results are robust for | < K < 1, and they re- 
fine the statement in Ref<21 that the two-channel Kondo 
ground state is stable for as long as K < 1. Our re- 
sults also provide the first theoretical realization of the 
quantum phase transition between 1CK and 2CK physics 
in Kondo impurity models. Further investigations via 
field-theoretical and Numerical Renormalization Group 
(NRG) 11 approaches are needed in order to clarify the 
critical properties, including the critical exponents and 
finite-temperature dynamics in crossover functions asso- 
ciated with the transition^. Our results could in princi- 
ple motivate the search for these critical properties near 
1CK-2CK quantum phase transition in future experi- 
ments on Kondo quantum dot coupled to 2D topological 
insulators. 
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Appendix A: The RG scaling equation in the 
weak-coupling regime 

In this Appendix, we provide some details on deriving 
the RG scaling equations of Eq. [3] for K —> 1~ from the 
bosonized Hamiltonian Eq. [3J Following Refs i 22 ' 23 , we 
decompose the boson fields <!>„ = 8 v ,<p v with v = s,a 
into the "fast" ) and "slow" components: 

$ v (t) = $<(r) + $>(r), 

\u n \<H' 

(Al) 
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with pf = p, + dp. The partition function can be decom- 
posed in the following form: 



D^<D^> e -So[<f < }-s [<f>]-s M l<s>>+>T < 



= Zn / ^$<e- s °[* < ] < e -^*t* < +* > ] >, 



(A2) 



where 



<A> f = J ^$ > e- So[$> U[$ > ]. (A3) 

The partition function Z^ can be re-expressed by expo- 
nentiating < ... >f in the integrand in terms of the effec- 
tive action S' e //[$ < ] = J drCx ( < & < ) with £ K being the 
Lagrangian of the Kondo model (see Eq. [2j, involving 
only the slow component of the fields with the following 
form via the cummulant expansion: 



5o[$ < ]-ln< e - s '-[* < +* > l >, 
< 5 mt [$< +$>] >/ 

5(<sg rt [* < + * > ] 



- <s mt [$ < 



$>] >}>) + 



(A4) 



The RG procedure is carried out by integrating out the 
fast modes of bosons and expressing the effective low- 
energy theory in the original form with the renormalized 
couplings. The following two-point correlation functions 
of boson fields prove to be useful in the RG analysis 22 : 



G(x,t) 



G(t) 



< $(x,r)$(0,0) >/ 
dk f du) 



2vr / 2tt 



— ikx iojr 

e e 



C( a T] - / ^ K ^' T ) fOT A » 1 
G(U,r) - < for / T i 



2tt 



(A5) 



where i^o is the Bessel function of the second kind. It is 
clear from Eq. IA5I that G(r) can be considered a short- 
ranged function of r. 

First, we focus on the first order cummulant < 
<Si„t[$ > + $ < ] >, which leads to the bare scaling dimen- 
sions of various Kondo couplings in Ref. 21 -. The renor- 
malization of the forward longitudinal term Jf term, 5 Jf , 
gives: 

Jf / dT<d x 8 s (0,T) > f 



= j{ \ d T [d x ef(o,T)+<d x e>(o,T) > f }. (A6) 



Since 9 S is an odd function in spin space, its average 
vanishes, < d x 9f(0 7 T) >/= 0, SJf = 0. This gives the 
first-order RG scaling equation: 



din p 







(A7) 



with the renormalized dimensionless coupling jf defined 
as: jf = poJ? where po — —77- is the density of states. 
The rescaling of backward longitudinal term Jf term 
leads to: 



/ drjf < sin( A /— a (O,r)) > f < sin(V2^0 Q (O,r)) > f 
A 



2K X 




J drJ|sin(^e<(O,T))sin(^rZ0<(O,r)) (A8) 

where Eq. IA5I and < e A >= e^ <A > are useo~,. Upon 
rescaling t, t — > T ^r^ we may define the new dimen- 
sionless renormalized coupling jf (p,) in terms of the bare 
coupling Jf (/1 = fio — 1) a s: 



j*(p) = Po p2+/ K 



(A9) 



we arrives at the RG scaling equation at the level of bare 
scaling dimension: 



d\n p 



= ( 



K+l/K 



-1)3 



2 • 



(A10) 



The first-order RG scaling equations for the remaining 
couplings are obtained similarly: 



dp? 
d\n p 

da? 

din p 



(K- 

,K 



l/K 



(All) 



with the renormalized dimensionless couplings defined 
in the text. Next, we consider the second order cum- 
mulant terms generated from — i(< Sf nt [^ K + $ > ]— < 
Smt[® < + < i )> ] > 2 >/)- In general, the second-order con- 
tributions to the renormalization of various couplings 
have the following form: 



dji 



xy 



d\n p 

djf 



■xy -z 

-aih fx 



d\np 

■ -XI 

h' 



■xy -z 

02J2 Jl) 



df 2 y 



dln/i 
djj 
dlrifi 



■xy -z 
-ClJi 32 



■xy -z 

C2J2 3i, 



o 7 -xy -xy 

= -2diJi" 3i , 



(A12) 



with dj, bi, Cj, and being the pre- factors to be deter- 
mined. 
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We first focus on the terms in J% v Ji which will con- 



tribute to the renormalization of j^ 1 
J2 J2 -S Z S+ I dr I dr' 



(Tra) 2 



/2tt 



x 2[<e-^^ ^cosU/^ Q (0,T)) 

V K 



x sin(W— a (O,r'))sin(V2^i?0 Q (O,r')) >/ 



< e 



-W2TrK<f> s (0,r) 



C0S(W— 0„(O,t)) >f 



We may use the following identities 22 to simplify Eq. lA16l 
V K 

= lim —d x < e'lVf^^^e-'^t^') 



x e -i(V2^K^<(O,T')) e iV^i r K0<(O,T') 



and 



(A17) 



'2tt 



x < sin(W— ^ a (O,r'))sin(V2^0 Q (O,r')) >/]. 

(A13) 

After averaging over the fast modes and rescaling t, r', 
we arrive at: 

J?" J? 



(7TO) 



x -l) e - t ^ F ^ < (°^cos(V2^0<(O,r')), 

(A14) 

In deriving the above equation, we have de- 
composed the terms sin(V 27T-K" (0 a (r' ) ) and 

cos(y / ^(6» a (r))sin(y / ^(6' a (T')) into the fast and 

the slow modes, and kept only the leading (more 
relevant) terms. In the limit of r, r' <C ~ a, we may 
get rid off one of the double time-integrals in the above 
equation by introducing a short-time cutoff tq ~ tt-. In 

the limit of if — > 1~, Eg. I A14I becomes: 



a x Vaz 

J2 J2 c+ 



.5+ / dr ^ e -iV^K0.<(O,r) COS (V2^X0<(O,T)). 

na J fj, 

(A15) 

Therefore, the pre-factor 02 in Eq. IA12I is found to be 
02 = 1. Similarly, we find the pre- factors c\ = b\ = 62 = 
di = 1 in Eq. |A"l2l 

Next, we consider a different type of renormalization 
involving Jf terms. We may focus on a typical term 
Ji V Jf , which renormalizes J^: 



2tt Jl y J{ 



dr / dT'S~S z 



K (ira) 2 
x [<M.(0,t) 

X e -»v / 2^f0 s (O,r') co^V^fi^o, r')) >/ 
- <9^ s (0,r)> / 



■ < ' ' > cos(\ / 2nK<f> a (0,T')) >/] 



-i\Z2TrK<p B (0,T') 

l2nJ* v Jf 



dr / dT'S~S z 



K (Tra) 2 
x [<9^>(0,r) 

X e -iV^0 8 (O,r') 008(72^^(0,/)) >/]• 



(A16) 



(A18) 



With the above relations, Eq. IA17I becomes: 



li m 1^ < e ^(^VTfl>(*,r)-^>(0,r')) > I 
»7->0 i?7 

X < e i >^*a(0,r') >/ 

v „-iv / 27Tff0?(O,T') o iv / 27rff0<(O,r') 

Ac- fc- . 

= 7(JflxG 9 > (jB)T) + 27r5 ai < 6f(x,T)<t>>(0,T') > / )| x=0 

x < e iV2^K0>(O,r') > ^ < el V27X0>(O,r') > ^ 

(A19) 

The leading logarithmic correction comes from the term 
d x < T )<Pf(0, t ') >/j which can be evaluated via 

the following relations^: 

i d(f) s 86 s 

v' F dr dx ' 

i d9 s d(p s 

v' F dr dx 



(A20) 



We have therefore 



-i d 



d x < 0>(x,t)4>(OS) >f = — 7r < 6>{x,t)6>{Q,t') > f 

Vp OT 



(A21) 



After collecting all the terms and performing re-scaling, 
Eq. IA16I becomes: 



1_ J J1R [trs-fax-^ 
v F ira J \x \i 



x e 



-i(v / 27T r ff0<(O,T')) 



cos 



(V2^<(oy)). 



(A22) 



Finally, the correction to j x contributed from J^J-f, 
reads: 

L '3l H (J, J 

x cos(v / 2lr¥</><(0,r)). (A23) 
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We therefore find a\ = 1. Similarly, we find c-i = 1. 
Combining the first and second order corrections to the 
renormalization of various Kondo couplings, Eq. [3] fol- 
lows. 



Appendix B: The 1-loop RG scaling equations near 
2CK fixed point via poor-man's scaling. 

In this Appendix, we derive the RG equations in Eq. 
fT3l from the effective Kondo Hamiltonian Eq. [6J in the 
weak coupling regime via poor-man's scaling as shown in 
Ref.— . Based on the scaling dimensions of various Kondo 
couplings in the strong couping 2CK regime, we take the 
logarithmic derivative of the proposed new dimensionless 
Kondo couplings j% v = poC2P e J 2 v ,ji = Pocfp e J* (see 
text) with respect to the cutoff energy p. 

First, we focus on the RG equation for j^ v : 



d In fi 



■ xy 

= % 



P Poc 2 



d In p 



The derivative of J 2 w.r.t. ln^t is given by: 



d,j*y 



d 



dhip In p 



r 1 du][j x 2 v j{ 

J Un 



xy J z p2xy,u(u)- l 



-Ul 



(Bl) 



(B2) 



Here, p 2x y,i z (^) is the effective electron density of states 
due to the additional phase correlations associated with 
the product of J% v and Jf terms in Eq. [51 Follow- 
ing Refi22, this is equivalent to replacing the free elec- 
tron Green's function of the effective leads: G a a =< 
ci a (t)c°(0) > by a "mixed" one: 



G' f w (t) 



\t)cf{Q) 



> 



x < e 



±n/47r(i-i)e a (t) 



>< S ± (0)S z (t) > 



(B3) 



(G^(t) can be defined similarly.) 

Therefore, p 2 xy,i z ( UJ ) = IT S CT Im ( G a( UJ )) reads^: 

P2xy,lz{u) = Po dEP 2 ± lz (E), 

Jo 

P 2 ±u{E) = ±-fdt< 6 2±lz (t) > e lEt (B4) 

where p = Im(G^ (uj)) is the constant density 

of states of the non-interacting leads, and < 2 ±i z (t) > 
has the following typical form: 



< OaxuW >=< e ± V 4 "(^- 1 )e«W >< S ± (0)S z (t) > . 

(B5) 

Since the exponential factors in 2 ±± z (t) are un-paired, 
it gives a trivial result: < e ±l vM¥^) e °(*) >= i. 



and hence it does not affect the renormalization of the 
couplings^. Nevertheless, < S ± (0)S z (t) > shows non- 
trivial correlations near—: 

< S ± (Q)S z (t) >« — (B6) 
with w c being a high-energy cutoff. We have therefore 



P2xy,l Z (u) = Po / dEP 2±lz (E) = p C 2 0J 
Jo 

1 



2e' 



(B7) 

with 5 2 = A(e') = ^sin(7re')r(l + e') with e' = ^ 
and r being the Gamma function. 
The integral in Eq. IB2I gives: 



J Mo 



-p 2 ctA(e')^/j^Jlln-^. 

Po 



(B8) 



Similarly, we find the RG scaling equation for jf is given 
by: 



dhip 



e'jl-plclA^^iJ?) 2 



(B9) 



where the effective density of states is used^: 
P2±,2±(u) = Po dEP 2 ±2±(E), 



P2±2±(E) = j dt< A(t) > e tEt , 



< A(t) > = < e -V 47r (i- 1 ))^( t ) e V 47r (-t- 1 )^( ) > 
1 



(BIO) 

We may determine the pre-factors c 2 , c z by the following 
identifications: = v / A(e)A(e / ), cf = A(e'). With the 
above results, we finally arrive at RG scaling equations 
shown in Eq |13l 



Appendix C: The 1-loop RG equations near 2CK 
fixed point via e-expansion technique. 

In this Appendix, we offer an alternative route to 
Eq. [T2] via e-expansion technique^— . 

We first derive the renormalization factors Zjx/z and 
Zf shown in Eq. [24j We focus on the 1-loop renormaliza- 
tion of the dimensionless bar couplings j% y = p^c^r p c J 2 V , 
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and jf = poclfi e Jf. Let us look at vertex renormaliza- 
tion of jf first— 



3i 



(Cl) 



where the effective density of states reads^: 

pj.(w) = po / dEP ± {E), 
Jo 



dt < A(t) > e %m , 



P ± (E) 



< A(t) > = < e -<V5f«.W e *VSFfla(0) > 



(iuj c t) 2 



(C2) 



From above, we find Pj_{E) = ci-E 2e_1 with ci = A(2e) 
and the constant A(e) being defined in Appendix B.. 
Therefore, we have 



2 c 



(C3) 



Plugging in these results into Eq. IC II and via the proper 
identification: c^" = y/c[ = -J 'A(2e), at the leading order 
in (jDVjf, X, reads: 



1 1 (jim 

2e 



Similarly, we can show that 



Zjs. = 1 



3i 



(C4) 



(C5) 



where the following relations are used: 
Po 



Pz{oj) 



dEP z {E) = pocau 



2c' 



2^ jo 

P Z {E) = J dt< B(t) > e lfct = c 2 E 

B{t) = <S ± (0)S z (t)>^-^-- r 

(iuj c t) e 



e'-l 



(C6) 



with C2 = A(e'), and cf = \fc\jci- 

Next, we provide derivation for Zf. Following Ref.— , 
the self energy at 1-loop order (see Fig. 5(a) of Ref—) 
leads to the following renormalization factor Zf for im- 
purity fermion: 



Z f = 1 



(C7) 

Plugging Eq. IC3l and Eq. lC6l into Eq. IC7l and expressing 
results in terms of the dimensionless couplings j% an d 
jf , we arrive at: 



Z\2 



z f = i + ^+ { j^r- (C8) 

With Zy, Zj-x , Zjz at hand, we now can reproduce the 
RG scaling equations Eq. [13] via the /3-function within 
the field-theoretical e-expansion approach : 



U2,0'-?1.0 



(C9) 



with the relations between the bare Kondo couplings jf , 
32 o ano - ^ ne renormalized ones jf , j 2 w being j| y = J% v = 



j 2 , and jf n = Jf = 
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